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Abstract. Surface growth models may give rise to instabilities with mound formation whose tipical linear 
size L increases in time (coarsening process). In one dimensional systems coarsening is generally driven by 
an attractive interaction between domain walls or kinks. This picture applies to growth models for which 
the largest surface slope remains constant in time (corresponding to model B of dynamics): coarsening 
is known to be logarithmic in the absence of noise {L{t) ~ \nt) and to follow a power law {L{t) ~ t^^^) 
when noise is present. If surface slope increases indefinitely, the deterministic equation looks like a modified 
Cahn-Hilliard equation: here we study the late stages of coarsening through a linear stability analysis of the 
stationary periodic configurations and through a direct numerical integration. Analytical and numerical 
results agree with regard to the conclusion that steepening of mounds makes deterministic coarsening 
faster: if a is the exponent describing the steepening of the maximal slope M of mounds {M" ~ L) we 
find that L(t) t": n is equal to j for 1 < a < 2 and it decreases from j to ^ for a > 2, according to 
n = a/{5a — 2). On the other side, the numerical solution of the corresponding stochastic equation clearly 
shows that in the presence of shot noise steepening of mounds makes coarsening slower than in model B: 
L{t) ~ t^/'*, irrespectively of a. Finally, the presence of a symmetry breaking term is shown not to modify 
the coarsening law of model a = 1, both in the absence and in the presence of noise. 

PACS. 68. Surfaces and interfaces - 81.10.Aa Theory and models of film growth - 02.30.Jr Partial 
differential equations 



1 Introduction 



In real systems surface growth occurs on two-dimensional 
(2d) substrates and therefore its modelization in one di- 
mension (Id) is in general an oversimplification, mainly 
justified by the possibility to have a deeper understanding 
of the dynamical evolution of the system. In some cases the 
surface indeed maintains a Id profile, for example when 
the relaxation of a grooved surface is studied H]. This is 
not true when the surface undergoes kinetic roughening 
or a growth instability followed by a phase separation: in 
both cases, noise makes the resulting morphology 2d even 
if the initial one is Id. The dynamical evolution of a vic- 
inal surface is another possible application of Id models. 
In fact, if atomic steps move in phase step motion is 
described by a one-dimensional growth equation. 

In this paper we are interested in a growing surface 
that undergoes a kinetic instability. The origin of such an 
instability is related to an additional barrier that diffus- 
ing adatoms must overcome to descend step edges (the so- 
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called Ehrlich-Schwoebel (ES) barrier |§,^). Even in the 
presence of stabilizing mechanisms, at sufficiently large 
scales the ES effect can be the dominant one and the flat 
surface becomes unstable against small deformations. In 
a first linear regime, a structure with a well-defined wave- 
length emerges and its amplitude increases in time. At a 
later stage, nonlinearities come into play and the mound 
structure typically -but not unavoidably coarsens. 

This simplified picture reminds spinodal decomposi- 
tion and coarsening that take place during phase separa- 
tion. Surface growth instability and phase separation do 
have strong similarities and they can indeed be equiva- 
lent in Id |q]. However, they are definitely different in two 
dimensions 

Generally speaking, phase separation has different prop- 
erties in one and two dimensions. We mention here two of 
them 10]: (i) In 2d, coarsening is driven by the tension of 
domain walls while in Id it is due to interaction between 
walls, (ii) Noise is generally irrelevant in 2d, while it mod- 
ifies the coarsening law in Id (except in the presence of 
long-range interactions) . 

The class of models that we study in one dimension 
is of interest in two respects. First, since the slope is as- 
sumed to increase indefinitely, it is not possible to speak 
of domain walls between different phases. Second, being 
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the interaction long ranged, noise may happen to be irrel- 
evant for the coarsening law in Id: this occurs for some of 
the models studied here. 

In the next Section we give a short introduction to 
the continuum equations that are encountered in one- 
dimensional models of conserved surface growth. A more 
detailed analysis can be found in recent review papers 
^ . In Section || we recall a few theoretical approaches to 
one-dimensional coarsening and in Section ^ we apply the 
linear stability analysis to our class of deterministic a- 
models. Its results are confirmed by numerics. The prob- 
lem of coarsening in the presence of noise is addressed 
directly via stochastic numerical integration in Section ^ 
In Section ^ we examine the effect of a symmetry breaking 
term, while the discussion of the results and our conclu- 
sions are presented in the final Section. 



2 Models for unstable surface growth 

In the Introduction wc have already mentioned the ES 
barriers as a source of the kinetic instability. In fact, it is 
well-known |^ that an asymmetry in the sticking coeffi- 
cients of an adatom to a step produces a slope-dependent 
current JesCti); h{x,t) and m = dxh are the local height 
and slope of the surface. It is more convenient to use the 
variable z{x,t) = h{x,t) — h, where h — Fgi is the aver- 
age height, instead of h{x, t) {F^ being the intensity of the 
flux). In this way the evolution equation dth = Fq — dxj 
simply writes dtz — —dxj while the definition of the slope 
m = dxh = dxZ is unaffected. 

For symmetry reasons, Jes is an odd function of m. 
Therefore we expect that Jes — t^Tn at small slopes, as con- 
firmed by a more rigorous analysis Q . The asymmetry in 
the sticking coefficients is generally due to the additional 
energy barrier hindering the adatom from descending a 
step. This gives rise to an up-hill current Jes — with a 
positive coefficient v. This implies that Jes has a destabi- 
lizing character, as easily revealed by the solution of the 
linear differential equation dtZ = —OxJes — —i^dxz{x,t) : 
z{x,t) = ^0 cos(gx)e'^''*, with ujq = vq^ . 

In a continuum model, where the lattice constant a 
goes to zero, the following expression for Jes can be used: 



model 1 . 



(1) 



If a discrete model with square symmetry is used, we 
expect Jes to vanish for m = toq = 1/a • We can therefore 
define the model: 



Jes — vm{\ — rn^ /m^) model 



(2) 



However, other mechanisms can produce a slope-de- 
pendent current: short-range step- adatom interaction 
pl3| , post-deposition transient mobility and downward fun- 
neling |ll| . The first mechanism may be either stabilizing 
or not, depending on the sign of the step-adatom interac- 
tion while the other (non-thermal) mechanisms are typi- 
cally stabilizing, i.e. they contribute with a negative term 



—v'm to Jes- Because of that, either the ES current ac- 
quires a zero at finite slope (from model 1 we pass to 
model 0) or there is a change in the value mo at which Jes 
vanishes. 

Finally, in a previous paper we have generalized 
Eq. (|^) to a class of models (a-models) characterized by 
different asymptotic behaviors for m — > oo: 



Jes 



(1 -t- ^2 „i2) 



2\a 



model a (> 1) . 



(3) 



In the Introduction we also mentioned possible stabi- 
lizing mechanisms. The simplest expression for a stabiliz- 
ing current is the so-called MuUins term that in its lin- 
earized form reads = Kd^m. Its origin may be either 
thermodynamic (relaxation through surface diffusion [p^ ) 
or kinetic (flu ctuations in the nucleation process of new 
islands ) • A stabilizing current gives a negative con- 

tribution to LOq. Starting from the equation dtZ = —dxij-Es + 
Jm) it is easily found that cug = lyq^ — Kq^. A flat surface 
is thereby linearly unstable against fluctuations of wave- 
length larger than Ac = i-KyjK/v. 

Both currents Jes and Jm change sign if a; ^ —x or 
z — !■ — z: the growth process cannot break the former sym- 
metry but it does break the latter. It has been shown 
|l^ that such symmetry breaking term (intrinsically non- 
linear) has the form Jsb = dxA^m^), where A ~ rn? at 
small slopes and A ^ —1/m? at large ones. Its presence 
is strictly related to the breaking of the detailed balance 
principle | [T^ and therefore to the non-equilibrium char- 
acter of the growth process. 

It has been proven jl^ that Jsb does not modify the 
coarsening law of model 0: the effects of Jsb on model 1 
will be considered in Section ||. 

We conclude this Section by writing down explicitly 
the class of growth equations that are analyzed in the 
paper: 

dtz{x, t) = -dxj{x, t) + ■q{x, t) (4) 



3{x,t) 



dim 



(1 

m(l 



model a 
m?) model 



(5) 



(r/(a;, t)) = (?7(x, t)?7(x', t')> = Fa5{x ~ x')5{t - t') (6) 

Eq. (^ is the evolution equation of the local height 
z{x,t) for a conserved growth process in the presence of 
the shot noise r]{x, t); Eqs. (||) give the currents for the a- 
models and model 0, once that x, t, z have been rescaled 
in order to get an adimensional equation; Eq. (^ gives 
the spec tral properties of the noise, whose strength Fq = 
Ft^aP / y/i>K is the only parameter appearing in the prob- 
lem: _Fo is the intensity of the flux, a is the lattice constant, 
£ is the diffusione length £d for a-models (Eqs. (0,||)) or 
the inverse of the constant slope mo for model (Eq. (0)), 
V and K are the prefactors of Jes and , respectively. 



3 Theoretical approaches to coarsening 

In this Section we review two theoretical methods that 
have been used to study the coarsening process in model 0. 
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The first method uses the property that the current Jes 
vanishes at a finite slope uiq . Although it can not be used 
for a-models, for completeness and clarity we discuss it in 
Section 3.1. The second method consists in a linear stabil- 



ity analysis of the stationary configurations. I t is a very 
general method and it is discussed in Section In the 
context of p hase separation processes it was introduced by 
Langer i20| to study model and its solution is proposed 
in Section [4.l| . Its application to a-models is carried out 
in Section \L2. 



3.1 Kink dynamics 

Model and a-models have the same linear behaviour, but 
they strongly differ in the nonlinear regime: in model the 
slope increases up to the maximal value tuq = ±1 while 
in a-models it grows indefinitely. 

For model 0, the surface profile corresponding to a con- 
stant value Too = 1 (i.e. to a vicinal surface with unitary 
slope) is stable but it can not be attained starting 
from a flat surface because the average slope must remain 
constant. Still, we can consider the stationary configura- 
tion TO+(a;) {j{m+{x)) = 0) corresponding to a limiting 
slope ±TOo for x ^oo. That profile is called 'mound' 
in the 2;-space and 'kink' or 'domain wall' in the space of 
the order parameter m, where it has the form to+(x) = 
tanh(2:/V2). 

During the coarsening process of model a typical 
surface profile is just an alternating sequence of kinks 
(to = TO_|_(a;)) and antikinks (to = m^(x) = — to+(x)), 
whose average distance L{t) increases with time because 
of the annihilation process between pairs of neighbour- 
ing kink-antikink. In the absence of noise the dynamics 
of kinks is governed by their attractive interaction that 
decays exponentially with the distance, since |TO-|-(a;)| ~ 
1 — 2exp{—^/2\x\) for |a;| ^ 1. Such a weak interaction 
determines a very slow coarsening: L{t) ^ Ini p9|,p0|. 

In the presence of noise the picture is different because 
of the induced fluctuations on the kink positions. If there 
were no constraint induced by the conservation of the or- 
der parameter, kinks would simply perform a random walk 
and therefore would travel a distance L in a typical time 
L^, giving a coarsening exponent n — 1/2 {L{t) ^ \/t). 
In a growth process, where a conservation law does exist, 
kink trajectories are not independent and the coarsening 
slows down: n = 1/3. This exponent is more easily under- 
stood in a spin picture |21|, where the conservation of the 
order parameter (the magnetization) implies that the sys- 
tem evolves through spin-exchange processes (Kawasaki 
dynamics). 

If we now turn to a-models, we recognize that the kink 
picture is not applicable because the unstable current Jes 
vanishes for infinite slope only. 




Fig. 1. Profiles of the different potentials Ves- From top to 
bottom: Model 1 (dotted line), Eq. (|VModel (dashed line), 
Eq. (1); Model a = 3 (full line), Eq. (|^. 



dtz{x,t) = 0, i.e. the current must be a constant: j = c. 
The net current c is related to the average slope of the 
surface: we are interested in a high symmetry surface and 
therefore we set = 0. 

The equation j = is formally equal to the Newton's 
equation for a particle of unitary mass, where the slope 
TO plays the role of the particle position and x the role of 
time: 

dlm{x) + j^s{m) ^ . (7) 

The fictitious particle feels the force — :/es(to), i.e. it 
moves in the potential Ves(™) = / dsj-Es{s). Different 
models give qualitatively different potentials (see Fig. |l|): 



VEs(m) = toV2- to'*/4 
Fes(to) = (l/2)ln(l + TO2) 



model 
model 1 



(8) 
(9) 



1/es(to) = -(1 + TO^)^-"/[2(a - 1)] model a > 1 (10) 

Stationary configurations therefore correspond to the 
periodic oscillations of the particle around the minimum 
of Ves in to = 0. We label the stationary configurations 
TO2l(x) through their period 2L. What about the limit 
L ^ oo ? For model it corresponds to the kink-solution: 
"^■00(2;) = m^(x) — *■ ±1 when x — *■ ±00. For model 1, the 
energy of the particle diverges when the period L increases 
and the limiting configuration m^{x) does not exist. It 
does exist for a > 1 and it corresponds to the well defined 
problem of a particle of zero energy moving in the poten- 
tial ([l0|): it starts at to = —00 and arrives at to = -l-oo 
after an infinite time. 

We start by considering small deviations ij) from the 
periodic profile: m{x,t) = m2-L{x)+%l)(x,t). Since Jes(w2l + 
'/') = Jes(™2l) + jEs('^2L)V' + C'(V'^): we obtain that the 
linearized evolution equation for iplx^t) is 



3.2 Linear stability analysis of stationary configurations 



(11) 



We now discuss the linear stability analysis of the station- 
ary configurations. They are determined by the condition 



and therefore the time dependence of -0 is: ')p{x, t) = (j){x)- 
exp(— et). Replacing the expression of '0 in terms of (j) into 
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Eq. (fi^j, we find that the stabiHty is determined by the 
spectrum of the following operator: 



(12) 



where = —d^ and H = —d'^ + Ui^{x) is a single-particle 
Schroedinger operator corresponding to the periodic po- 
tential Ui^{x) = — j^g(m2L(a;)), of period L: 



U^{x + L) = -j'^^{m,^{x + L)) 

= -Jes(-"^2l(2;)) = U^{x) 



(13) 



In one dimension coarsening is due to the unstable 
character of the periodic stationary configurations, i.e. to 
the existence of negative eigenvalues in the energy spec- 
trum . Because of the periodic character of the opera- 
tor DxH, eigenvalues are grouped into bands. 

Our evolution equation for m(x, t) is dtm = D^j where 
the current j (see Eq. (|^)) can be derived from a pseudo 
free energy T: 



J 



5T_ 

5m 



— T 



dx{\{d^m.f ^V^^{m)\ (14) 



For model the potential —V^^(rri) has the standard 
double well shape and m(x, t) evolves accordingly to the 
Cahn-Hilliard equation; in the presence of conserved noise 
we obtain the so-called model B of dynamics |2^. If 
is replaced by the identity operator, the order parame- 
ter is no longer conserved and its evolution equation is 
dtm = j, which is equivalent for model to the time de- 
pendent Ginzburg Landau equation, or -in the presence of 
nonconserved noise- to model A of dynamics . We use 
the notations e, 4> for the spectrum of the operator D^H 
and e, for the Hamiltonian operator H. 

Let us start with few general statements on the lowest 
part of the spectrum. Translational invariance implies that 
e = is always an eigenvalue of the operator H (and 
therefore e = is an eigenvalue of D^H as well) and it 
corresponds to the eigenfunction 4>(x) = 4)(x) = m'^^{x). 
To prove it let us use the definition of m^i^{x) as solution 
of the differential equation (R) and take its derivative: 



TT^'lix) + JEs(n^2L)m;ja;) = 
Since Ui^{x) = — JEs('^2L(a;)) we just have 



(15) 



- a^(m;,(x)) + U^{x){m',,{x)) ee Hm',,ix) = . (16) 

Therefore the operators H and DxH have a zero en- 
ergy eigenvalue and the corresponding eigenfunction m'^^ {x) 
has period X — 2L (i.e. it corresponds to the wavevector 
q = tt/L in the Bloch representation). 

We can now recognize the importance of the limit 
m^{x). If moo (x) exists, for L — > oo the periodic potential 
J7l (x) becomes a single well potential Ur^ (x) . In this limit 
m2i^{x) is a monotonic function rn_|_(a;) (the kink solution, 
for model 0). Therefore 4ii{x) = (/>i(a;) = m'_^{x) has no 




Fig. 2. The single well potential Urx,{x) for model (dashed 
line) and for model a — 3 (full line). At large x the former 
approaches 2, while the latter vanishes as U^{x) = a/x"^ , with 
a — 10/9 (dashed-dotted line, see formula (123)). 



node and it represents the ground state for the single well 
problem. For finite L we have a periodic potential and the 
energy level ei = ei = gives rise to the lowest band of 
the spectrum. The ground state €gs{L) of the operator H 
corresponds to g = and must therefore have a negative 
energy, implying that H has negative eigenvalues. The re- 
lation Cgs = ^{l — 0) is no longer valid for the conserved 
model, but Langer has shown that negative eigen- 
values e of if may be put in correspondance to negative 



eigenvalues e of D^H (see Section 4.1). 

Since an unstable mode increases exponentially with 
the factor exp(|e|t), the knowledge of the L dependence of 
the ground state energy allows to find the deterministic 
coarsening law L{t) via the relations |eGs(i)| \/t (non- 
conserved model) and |eGs(i)| ^ 1/t (conserved model). 



4 Deterministic coarsening 

In Fig. M we plot the single well potentials U^{x) for 
model (dashed line) and for model a = 3 (full Hne), the 
latter being representative of all the class of a-models. 
Since, as already explained, ei = is the ground state 
energy for the single well, completely different behaviours 
are expected in the two cases. 

For model (see Section 4.1), U^{x) approaches 2 at 



large x: therefore the wavefunction 4>i{x) decays exponen- 
tially at large distances and energy shifts due to the tun- 
neling between wells at finite distance L are expected to 
be exponentially small. This property is the counterpart 
of the exponentially vanishing interaction between kinks, 
discussed in Section 

For a-models (see Section [i!^), ?7«,(|a;| = oo) = = ei 
and the considerations developed for model do not apply. 
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(a) 



1 - 




-3 -2 -1 



Fig. 3. The true periodic potential Ui,{x) (empty circles) is 
compared with the potential Ul{x) (full line), obtained as a 
superposition of the single well potentials Uaa{x — nL) (see 
Eq. (USl)). We make this comparison for model (a), where 
L ~ 6.8, and for model a —3 (b), where L ~ 22.8. In both 
cases we display one period only. Ui,{x) is obtained numerically 
and Ul (x) is obtained analitically. The superposition principle 
works well in (a) even for small L, while the exact potential 
for model 3 (see b) is not reproduced by the sum of single 
well potentials. In order to make more evident the difference 
between the two potentials, in (b) we focus the plot on the 
region of positive Ui,{x). We remind that Ui,{0) = — 1 for any 
a. 



4.1 Model 

For model the kink-solution is m^{x) — tanli(x/\/^, 
the single well potential is (see the dashed line in Fig. 0) 

U^ix) = -l + 3ml^{x) = -l + 3tanh^(x/V2) (17) 

and the corresponding ground state wavefunction is (f)i (x) = 
'«L(^) = 7fSec2(^). 

For finite L, Ui^{x) is a collection of wells centred at 
points X = 0, ±I/,±2L, ... . It is interesting to compare 
Ui,{x) with the periodic potential U*{x) obtained as a su- 
perposition of the single well potentials U^{x — nL): 



U:{x)^UUx) + J2[UU^ 

n#0 



iL)-C/^(w)]. (18) 



For model 0, U^{oo) = 2 and the summation can in- 
deed be limited to the two terms n = ±1 because the quan- 
tity in square brackets is exponentially small. In Fig. Ha we 
show that U*{x) (full line) is an excellent approximation 
of C/l (x) (circles) . 

Langer performed a tight-binding-approximation to de- 
termine the lowest band of the energy spectrum arising 
from the ground state ei = of the single well. The re- 
sult l2Cll is 



e{q) ~ —(1 + cosqL) exp(— 2L) , 



confirming that e{q = ^r/i) = 0, Cqs = £(? = 0) and 
that ecs decays exponentially with the distance L between 



wells. The relation jecsl ^ l/t implies that coarsening 
is logarithmically slow: L{t) ^ \nt. This result is valid for 
the conserved model as well. Langer proved it by means 
of the variational condition 120(1 : 



where 



l) 



dx(l)l{x)4>i{x) 



(20) 



(21) 



is the scalar product between the ground state function 
01 (cc) of the Hamiltonian operator and 4>i is defined by 
the relation Dx(f>i = 4>i. 

Because of the conservation law, the lowest energy 
band for the operator D^H reads [M: 



( . . 2 r exp(-2£) 
€{q) ~ - sm qL 



(22) 



The factor L at denominator is irrelevant and the coars- 
ening law L{t) ^ ln< is still valid. Eq. ( p2|) also confirms 
that e{q = tt/L) = and that e(g = 0) is no more the 
ground state. 



4.2 a-models 

The solution of model by Langer is of great interest 
because it gives approximate expressions for the lowest 
energy band, both for the nonconserved ( p^ ) and con- 
served ( |2^ ) models. The basis of his treatment is the tight- 
binding-approximation: AU{x) = U*{x) — Ura{x) is taken 
as a small perturbation of U^{x). 

In the case of a-models, a simpler strategy can be fol- 
lowed |l2j by replacing the periodic potential C/l (a;) with 
a double well U2{x), obtained by joining rather than su- 
perposing (x) and (x — L) ||2^ . Using this approxi- 
mation, our analytical results for the coarsening exponent 
n{a) agre well with the estimate obtained from numerical 
integration of the growth equations |E2| (see the results for 
the conserved model reported in Fig~P). In the following 
we report the general lines of this theoretical approach, 
and we add a direct numerical confirmation (see Fig. ^). 

The single well potential U^{x) for the a- models (see 
Fig. 1^) decays to zero at large x from positive values. From 
the relation U^{x) — —j'j^g{m^{x)) it is simple to derive 
that for TO ^ 1: 



{2a - 1) 
TO^"(x) 



(23) 



The integration of Newton's equation (^ for zero en- 
ergy gives the asymptotic result {x ^ 1): 



(19) which implies 



mZix) 



{2a- l){a- 1) 1 



a 



(24) 



(25) 
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Table 1. The deterministic coarsening exponent n, for the 
nonconserved (D^ = 1) and the conserved {Dx = —d^) models. 
For Q = 2 there are logarithmic corrections (see jl2|). 



1 < Q < 2 a > 2 



nonconserved 
conserved 




Thus the single well potential decays as the inverse of 
a square law whatever is a, but with a prefactor a that 
increases between a = for a = 1 and a — 2, for a = oo. 
For a = 3, a — ^ and the corresponding function ^ is 
reported in Fig. g as the dashed-dotted line, showing the 
comparison between the asymptotic expansion (|25| ) and 
the exact expression p^ . 

The solution (piix) of the Schroedinger equation for the 
single well (let us remind that the ground state has zero 
energy) therefore decays at large a; as a power law |1T2|: 



(j}i{x) 



^, with P = {1 — a ^). So, the ground state 



wavefunction is a bound state for > ^ i.e. a > 2 only. 

In Ref. jlj] we used the Landau and Lifshitz approach 
p?! for the double-well problem and we extended it to 
take into account the possibility that the ground state 
wavefunction for the single well is not a bound state. This 
happens for 1 < a < 2. The main point is that even 
if (j)i{x) is not a bound state, the ground state (f>2{x) of 
the double well is bound, because its energy £2 is now 
strictly negative and so lower than the asymptotic value 
C/2(oo) =0. 



Following the above approach we found |12 that 



7 = 



2 



a <2 
a>2 



(26) 



Once 7 is known, the coarsening exponent n is just 
given by n = 1/7. In Table |l| we summarize the results 
found in Ref. |]l2| for the nonconserved and conserved 
models. For a — 2 there are logarithmic corrections: L ~ 
(t/lni)". 

The reason why it has not been possible to treat the 
periodic potential Ui^{x) in a more rigorous way is given 
in Fig. pb, where it is clearly shown that the superpo- 
sition principle does not work for a-models: accordingly, 
the potential Ui^{x) can not be approximated as the sum 
of single-well potentials. Therefore, the application of the 
tight-binding-approximation is not straightforward becau- 
se we do not know the explicit expression of the perturba- 
tion AU{x) = Ui,{x) — Uoo{x). The disagreement between 
Ui^{x) and Ul{x) might appear at first sight to be unim- 
portant, but if we neglect their difference we obtain wrong 
results (as we have verified). 

However, it is possible to check in a direct way the ac- 
curacy of our theory: we solve numerically the quantum 
mechanical problem to determine the ground state energy 
eGs(i) of the full periodic potential Ui,{x) (see Fig. ^ and 
footnote psi ). The numerical results for the exponent 7 
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Fig. 4. Nonconserved models: on a log-log scale (main) we re- 
port the absolute value of the ground state energy ecsiL) for 
the periodic potential U-l{x). Circles refer to model 3 and dia- 
monds to model 10. In the inset we report the local derivatives 
d[log(— eGs)]/rf[log(-£/)] and the asymptotic values of 7 that are 
calculated analytically (see Eq. (pq)). 



agree fairly well, at large L, with the theoretical predic- 
tions 7(3) = 2.3 and 7(10) = 2.8. This confirms that our 
theoretical approach to calculate e^siL) is indeed correct. 



4.3 Numerical analysis 

Let us now discuss our numerical results for the determin- 
istic a-models. We have numerically integrated the equa- 
tion of motion (^ with 77 = 0, starting from an initial 
profile z{x,0) — r^, where r^ is a random variable with a 
flat distribution in the interval [—0.1,0.1]. 

Tipically, we have followed the dynamical evolution 
for a total time ^max ~ 400, 000 - 1, 600, 000, for a chain 
length Lc = 1024, with a spatial resolution 5x = 0.25 and 
an integration time step r = 0.05. A few tests have been 
also performed with a smaller time step r = 0.025 and 
with longer chains (Lc = 2048 — 4096), obtaining consis- 
tent results. The adopted integration scheme is a time- 
splittin g pse udo-spectral code: more details are reported 
in App. AI] . 

On top of Fig. H we display a portion of the surface 
profile z{x) for the model a = 3. It appears as made up of 
constant slope regions separated by domain walls. How- 
ever, the slope profile m(x) reported in the centre of Fig. |^ 
does not corroborate this picture: no region of constant 
slope is clearly visible and the maxima or minima have 
appreciably different values. The same remark is applica- 
ble at later times. On the bottom of Fig. |^ we also display 
the potential U{x) = —j^^{m{x)) that enters in the ana- 
lytical solution of the problem (see the previous section). 
It is therefore reassuring that U{x) looks indeed as a regu- 
lar sequence of the single well potentials depicted in Fig. |^ 
as a full line, because it confirms that the surface profile 
keeps close to a stationary configuration. 

The next step is to evaluate the characteristic length 
L{t), corresponding to the average distance between wells. 
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Fig. 5. Model a = 3, late stages of coarsening. Top: the sur- 
face profile z{x) as obtained via numerical integration. Cen- 
tre: the slope profile m{x) = z'{x). Bottom: the potential 
U{x) = — iEs("^)- We display only a piece of the total spa- 
tial domain. 
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Fig. 6. Deterministic coarsening: L vs t in log-log scale for 
a — 1.5 (o), Q = 3 (□) and a = 10 (A). For presentation 
purpose the data for a = 3 and 10 are shifted by a constant. 
The lines indicate the best fit to the data for t > 10, 000 and 

the slopes are equal to 0.250 (— ), 0.233 ( ) and 0.208 (• ■ ■), 

in agreement with the theoretical prediction (see Table qI). 



We define the wavevector Pc{t) via the relation: 



Pc{t) 



(27) 



where S{p,t) = |z(p, t)p is the power spectrum associ- 
ated to the field z at time t {z being its spatial Fourier 
transform) and the sum is restricted to the wavevectors p 
for which S{p) > S ■ Sm, Sm being the maximum value of 
the spectrum and 6 some threshold (typical values are 
5 ~ 0.1 — 0.2). The characteristic length is then eval- 
uated as L(t) — 2TT/pc(t) and the coarsening exponent 
n{a) has been obtained by considering the scaling be- 
haviour of L{t) ^ t" in a time interval 10,000 < t < 
400,000- 1,600,000. 

As an independent check we have also determined L 
from the normalized spatial correlation function of the 
surface profile 



C(r,t) 



< z{x + r, t)z{x, t) > ~ < z{x, t) >^ 
< z{x,t)'^ > - < z{x,t) >2 



(28) 



where the spatial average < • > is performed along the 
chain. Defining L through the relation 



C(L,t) =C(0,i)/2 



(29) 



our results are in agreement with the previous ones ob- 
tained from the power spectrum in all the considered cases. 

The numerically estimated length L{t) is reported in 
Fig. for a — 1.5,3,10. The coarsening exponents are 
n(1.5) = 0.250 ±0.003, n(3) = 0.233 ±0.005 and n(10) = 
0.208 ± 0.002. These results are consistent with the theo- 
retical estimates for n{a), as summarized in Fig. |[ Their 
critical discussion is deferred to the final section. 



5 Coarsening with conservative noise 

The stochastic equation (^ has been integrated up to a 
time iMAx ^ 800,000 with = 1024. The results for 
L{t) and C(r, t) have been obtained by averaging over Nc 
different initial conditions with different noise realizations 
(tipically A^c = 10). The integration scheme employed in 
the noisy case is different from the one adopted for th e de- 
terministic case and it is described in detail in App. |A.2| . 
We have used a noise strength corresponding to a value 
of Fq (see Eq. (^)) equal to 0.05. This is a physically rea- 
sonable value because for large ES barriers §, -Fb ~ 
and £u is typically of the order of a few dozens of lattice 
constants. 

For the noisy case we have defined L only in term of the 
average correlation function C(r, t), where the bar means 
that average is now performed at each time not only along 
the chain (see definition (|28|)) but also over different noise 
realizations. A sufficiently good scaling is obtained in the 
time interval 40, 000 <t< 800, 000 for all the considered 
values of a (= 1, 2, 10). 

As a benchmark to verify the validity both of our in- 
tegration scheme and of our procedure to estimate n{a), 
we have analyzed model 0. In this case the coarsening ex- 
ponent is known ||l^,^ to be n = 1/3. A good agreement 
between our numerical data and the theoretical prediction 
is found for t > 40, 000, as shown in Fig. 0. 

For the two values of a, a = 1 and a = 10, we find 
n = 1/4. We conclude that, in the presence of shot noise, 
the coarsening exponent is independent of a and equal to 
1/4. Fig. 1^ also suggests the possible existence for model 
of an intermediate regime (with L ~ 5 — 20) where an 
effective exponent n = 4- is found. 
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Fig. 7. Coarsening in the presence of shot noise: L vst m log- 
log scale for model (O) and for a = 1 (A) and a = 10 (o). 
Dashed lines have slope 1/4 and full line has slope 1/3. 



6 Effects of a symmetry breaking term 

In Ref. [0 one of us studied the effect of symmetry break- 
ing on model 0. Since in that model the slope keeps finite 
with a maximal value equal to one, the detailed expression 
of the function A{m^) {jsB = d^A) is not relevant. It was 
therefore chosen the simplest form, the one valid at small 
slopes: A(m^) = A*m^. 

On the contrary, for a-models the slope can diverge so 
that the exact expression of Jsb should be used |14|: 



1 



1 



(30) 



We limited ourselves to the physically relevant case 
a — 1. We have therefore integrated the following differ- 
ential equation: 



dtz 



-da: 



dim 



1 -|- 



1 



for values of A* varying between 0.1 and 1. 

Our results (see Fig. ^) suggest that jsB is irrelevant 
for the coarsening law, both for 77 = (deterministic case) 
and for rj ^ (noisy case). 



7 Discussion and conclusions 




Fig. 8. Coarsening for the asymmetric model (Eq. with 
A* — 1), in the absence (empty circles) and in the presence 
(full triangles) of shot noise. Fits have been done for t > 10* 
and give n — 0.25 without noise (dashed line) and n — 0.24 
with noise (full line). Circles have been shifted by a constant. 
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Fig. 9. The coarsening exponent n as a function of a, for the 
deterministic models (full circles) and for the stochastic models 
(open diamonds). In the inset we enlarge the region of small 
a. Full line is the theoretical result in the absence of noise 
(Table ^ and the dashed line is our ansatz n = 1/4 for the 
noisy case. 



In Fig. ^ we have summarized our numerical and theoret- 
ical results for the coarsening exponent n (L ^ t"). In the 
absence of noise, our theory (full line) predicts that n — j 
for 1 < a < 2 and for larger values of a it decreases down 
to ^ {n — 1/(5— ^)). Numerical results (full circles) agree 
well with the full fine. 

We are aware of only one analytical paper treating 
our class of models (Ref. [|2^). The author uses scaling 
arguments to conclude that, in the absence of noise, n = 
1/4 irrespectively of a and for any dimension d of the 
substrate. In the following we give a drastically simplified 
version of the scaling arguments. If Z, L and M = Z/L are 



respectively the typical height, width and slope of mounds 
at time t, the evolution equation for z{x,t) implies Z/t ~ 
j /L. The current j is made up of the MuUins term, of order 
M I L"^ plus the ES current, whose asymptotic expression 
for large slope is 1/M^"~^. They vanish in the limit t ^ 00 
and must be of the same order in which entails the 
relation M" ^ i. If their sum j is supposed to be of 
the same order as well, the relation Z/t ~ M/L^ implies 
Lit) ^ t^/^, i.e. n — 1/4 for any a. 

One main drawback of the scaling argument is that 
the two terms appearing in j are of the same order but 
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their sum is smaller (i.e. of higher order in j) because of a 
compensation effect. This is a necessary condition for the 
stationary configurations to play a role in the coarsening 
process. 

In order to have a direct numerical check of our state- 
ment, we have evaluated the quantities {\dxjM\), {\dxjEs\) 
and {\dx{jM + jEs)\), as a function of time, where (• • •) 
means, as before, the spatial average. For all the consid- 
ered values of a (a = 1, 3, 10), the result is the same: the 
ratio (|i9a;jM|)/ {\dxjEs\) is equal to one, up to higher order 
terms, and {{OxUm + jEs)\) / {\dxjm\) vanishes. 

Our numerical results tell even more than that: in fact 
scaling arguments would suggest that n is strictly smaller 
than J if j is smaller than Jm and j^s, but for a < 2 we 
do find 71 = i even if IjI/IjmI ^ bl/ljEsI ^ 0. 

Model a — 1 had been previously studied numerically 
at short times also in Ref. and authors found a value 
n w 0.22, independent of the noise strength. 

Let us now discuss the results in the presence of noise. 
Fig. |9| presents with diamonds the numerical results for 
the stochastic integration of Eq. (^) . Our data refer to a = 
1, 2, 10 and provide a reasonably convincing evidence that 
in the presence of noise the coarsening exponent remains 
constant, n = \ (dashed line). The somewhat larger value 
for a — 2 may be due to unknown logarithmic corrections. 

Authors in Refs. |^,|2j use qualitative arguments to 
describe coarsening assisted by noise: they use a 'single 
mound' model and find the coarsening time by requiring 
that shot noise induces a height fluctuation of the same 
order of the mound height. In one dimension they find n = 
1/(3 + ^), where a is defined phenomenologically through 
the asymptotic relation M ~ L^/" between the typical (or 
the maximal) slope M and the width L of mounds. 

Their prediction for n(oo) seems to agree with the re- 
sult n = i for model 0. This is reasonable because in 
that model the slope is constant and therefore it can ef- 
fectively be equivalent to the model a = oo. Actually, if 
we take the limit a = oo in Eq. (^), it is straightforward 
to conclude that the current Jes vanishes and we obtain 
the linear equation: 



dtz{x, t) = —dxz{x, t) + rj{x, t) 



(32) 



the current Jes is irrelevant provided that the slope m 
diverges. 

We can now summarize our main results. Without 
noise, n = ^ for a < 2 and n — 1/(5 — ^) for a > 2. 
This result has been obtained analytically and it has been 
confirmed by extensive numerical calculations. It can not 
be deduced by simple scaling arguments. In the presence 
of noise our numerical data for a = 1,2, 10 suggest that 
n = J irrespectively of a. This guess agrees with the well 

known result ^ {t) ~ t^^^, valid for the linear model a = oo. 
So, steepening of mounds makes coarsening faster without 
noise and slower with noise. 

We believe that the surface profile can not be described 
as a sequence of mounds with a spatially constant slope 
AI that increases in time (see Fig. centre). This wrong 
assumption may be the reason why qualitative arguments 
to determine n do fail. 

We have also considered the possible effect of a sym- 
metry breaking term Jsb in the current: it is irrelevant for 
model 1, as already proved for model Q. 

We conclude the paper by mentioning a different model, 
whose coarsening properties bear some similarities with 
our a-models. It has been studied by Bray and Ruten- 
berg and it consists in the addition of a long-range 
attraction between kinks to model 0. If such interaction 
decays as a power law of the distance {l/L"^ with s > 
1) the deterministic coarsening exponent is found to be 
n = 1/(1 -I- s). In the presence of conservative noise, this 
appears to be relevant for s > 2 and in that case n = ^. 
Some analogies therefore exist with our class of a-models, 
because in both cases the coarsening exponent is a contin- 
uously varying function of a parameter (a or s), noise may 
be relevant (a, s > 2) or not (1 < a, s < 2), and -finally- 
the stochastic coarsening exponent is constant if noise is 
relevant. 
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Lepri and M. Moraldi. 



The basic question is whether our class of a-models 
tends -in some sense- to Eq. ( |32| ) with increasing a. In 
the absence of noise, the answer is surely negative: indeed, 
Eq. (|3^ ) does not admit stationary periodic solutions that, 
as discussed above, are crucial for deterministic coarsen- 
ing. 

On the other hand, if noise is present (77 ^ 0), Eq. ( ^ 
describes a process of kinetic roughening: the growing sur- 
face is characterised by a correlation length S.{t) ^ t^^^, 
where z is the dynamical critical exponent. In d = 1 it is 
well known Q that for the quartic linear equation (^^, 
z = 4. It is reasonable toguess that the stochastic a-model 
does converge to Eq. (p2) for a — *■ 00 and -in the same 
limit- n{a) ^ \ = \- 

The meaning of a constant value n = j for any a 
is simple: in the presence of noise the detailed form of 



A Integration algorithms 

Let us rewrite in an explicit way the evolution equation 
(|) for the field z{x,t): 



dtz{x, t) = —dxZ — dx 



dx 



{l + {dxZY)" 



+ V{x,t) , (33) 



where ri{x,t) indicates additive (5-correlated spatio-tempo- 
ral gaussian noise, i.e. 



< r]{x,t) >^ 
<?7(a;,i)r/(0,0) >^ Fo5{x)S{t) 



(34) 
(35) 
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A.l Deterministic equation 



A. 2 Stochastic equation 



Let us first neglect the noise term: in order to perform 
the numerical integration of ( |33| ) we consider a discrete 
spatial grid of resolution 6x and a discrete time evolu- 
tion with time step r. The discretized field is written as 
z{i,n), where the integer indices i and n are the spa- 
tial and temporal discrete variables, respectively. Periodic 
boundary conditions have been considered for the field: 
z{i,n) = z(i + I,n), where / is the number of sites of the 
grid {Lc = ISx). The algorithm adopted to integrate 
is a time-splitting pseudo-spectral code |Q. In particular, 
by following Eq.(^3|) has been rewritten as 



dtz{x,t) = {£+J\f)z{x,t) 



(36) 



where £ and Af are two operators defined in the follow- 



ing way: Cz 



-diz and Mz 



As 



usual for time splitting algorithms, the linear evolution, 
ruled by the operator £, is treated independently from 
the nonlinear one (associated to the operator Af) . A com- 
plete evolution over an integration time step r therefore 
corresponds to the two successive integration steps: 



t + t) = exp [£t]z(x, t) 



and 



(37) 
(38) 

(39) 



Eq.(^9|) can be easily solved in the Fourier space and 
the equation of motion for the spatial Fourier transform 
of the field z{p, t) is 



z{x, t + r)— exp [N't\z*{x, t + t) 

where z*{x,t) is a dummy field. 

Let us firstly consider the linear part, 

dtz{x, t) = Lz{x, t) . 



4~ 

-p z 



dtz{p, t) 

The time evolution for z is simply given by 
z*{p,t + T) = exp[-p*T]2(p,t) 



(40) 



(41) 



Therefore, in order to integrate Eq.(|39|), the field should be 
Fourier transformed in space {T), then multiplied by the 
propagator reported in Eq.([fl|) and the outcome of such 
operation should be finally inverse-Fourier transformed 
(^-1): 



z*{x,t + t) = ^ exp [—p'^t\J^z{x^ t) 



(42) 



The integration of the nonlinear part has been per- 
formed by emploing a second-order Adam-Basforth scheme 



z{x, t + r) = z*{x,t + t) + 



(43) 



-[3Gizix,t))-G{zix,t^T))] 



where G{z{x,t)) ~ |— 9j; 



__d^_z_ 



In order to ob- 



tain a better precision, the spatial derivatives appearing 
in G{z{x,t)) have been evaluated in the Fourier space. 



Let us now consider the noisy problem: in this case the 
algorithm outlined here above does not guarantee a suffi- 
cient precision. Therefore, we have developed a more ac- 
curate integration scheme that consists as a first step 
in rewriting the equation of motion (|3^) in the Fourier 
space 

dtS{p, t) = -p^-zip, t) -t- Gp(z, t) + flip, t) , (44) 

where Gp{z,t) and f}{p^t) are the Fourier transforms of 
the nonlinear part and of the noise term appearing in 
Eq. (|33|). The amplitudes f]{p,t) of the noise components 
in the Fourier space are still gaussian 5-correlated stochas- 
tic variables with zero average and with a variance Vp = 
J-^/ 1 independent of p (white noise). The formal exact 
solution of is 



(45) 



t+r 



i{p,t)+ / dt'eP"^''-'\Gp{z,t')+f,{p,t')) 



z{p, t 



The problem is now to evaluate the two terms ap- 
pearing in the integral. The first term has been evaluated 
adopting a second order Adam-Bashfort scheme, i.e. 



t+r 



di'ef'(*'-*)G'p(z,i') = 
3Gp(z, t) - e-P'''Gp{z, t - r)] + 0{t^) 



(46) 



The treatment of the second term is more delicate, 
since it is a stochastic integral: we have chosen to evaluate 
it accordingly to the Ito's prescription | |37{ : 

dt' eP '^''-'^{p,t')^ (47) 



[Wp{t + T)-Wp{t)]^ AWp{t). 

Here Wp{t) and AWp{t) represent two Wiener pro- 
cesses: in particular we have {AW„) — and {{AWpY) — 
rVp. The complete solution of ( [f4|) can be written as 



z{p,t + r) = e P 
3r 



z{p,t) 

T 



(48) 



—e-p"^Gpiz, t) - -c-2p'^Gp(z, t-T)+ e-p"^AWp{t) . 

In order to obtain the solution in the real space it is suf- 
ficient to inverse-Fourier transform t + r). 

Due to the spatial and temporal discreteness of the 
integration scheme, the spatio-temporal noise term r/(a;, t) 
should be rewritten as ^f, where 7 is a random gaussian 
variable of zero average and variance 



V = 



Sx 



(49) 



with 



< l?lo >= Vd.oSn.o ■ (50) 
The spatio-temporal discrete gaussian noise, with zero 
average and standard deviation \/ has been numeri- 
cally generated by employing a Box-MuUer algorithm [Q . 
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